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We study the mutually coupled, strongly interacting bosonic and non-interacting fermionic, species 
of unequal masses in the regime were the retardation efltects are an important part of the physics. 
A cloud of neutral atoms experiences a synthetic magnetic field because of a vector potential that 
imposes a phase shift on the constituents. The magnetic field causes the oscillations of the magnitude 
and sign of the effective interaction between bosons from repulsive to attractive in contrast to the 
static case. We show that the dynamics for the gaseous Bose-Fermi mixtures when reaching the 
quantum-Hall regime becomes highly nontrivial. 
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The ultracold gases of atoms allow for the observa- 
tion and control of many-body quantum effects at macro- 
scopic scales. They can thus offer a fertile arena for 
exploration of condensed matter physics. The obvious 
limitations that come from neutrality of atoms preclude 
the observation of great variety of fundamental phenom- 
ena i.e. charge particle moving in magnetic field. How- 
ever, the analysis of the latter can be carried consider- 
ing rotating Bose-Einstein condensates [HE] trapped in 
lattice potential created by lasers. In a frame of refer- 
ence rotating about the z-axis with angular velocity 
the kinetic term in Hamiltonian is equivalent to that of 
a particle of charge Q experiencing a magnetic field B 
with QB ~ 2TOf2, where m is the mass of the particle 
This connection shows that the Coriolis force in the 
rotating frame plays the same role as the Lorentz force 
on a charged particle in an uniform magnetic field fS^, '3] . 
The above setting comes with limitations because a large 
magnetic fields / = mo?VL/i:h (angular momentum) are 
required to make possible the study of poorly explored 
bosonic states in case when / = p/q, {p and q is the 
ratio of atom number to the number of flux quanta re- 
spectively) is a rational number. Very recently Lin et al. 
[7] circumvented the problem by imprinting a quantum 
mechanical phase on the neutral atoms without rotating 
them at all and the quantum-Hall physics can be reached. 

Current experiments jHHH] on trapped mixtures of the 
atomic Bose-Fermi (BF) and Bose-Bose gases show that 
the presence of a relevant fraction of one modifles the 
quantum phase transition occurring in the other induc- 
ing a significant loss of coherence. These observations 
are supported by a theoretical description that includes 
the multiband virtual transitions [13], different masses of 
strongly interacting particles fT31 and numerical calcula- 
tions [15j. The density-density (DD) interaction between 
different species can be repulsive or attractive and is pro- 
duced by changes of one species density that induce a 
modulation of another. Therefore the dynamics underly- 
ing the phase transitions in the BF mixtures is produced 
by the feedback of the density perturbation and a shift 



of the inter-bosonic potential occurs, that changes the 
original interaction between them [16J providing various 
novel phases |17) . 

In the present paper, motivated by recent experiments 
done by Lin et al. [71 , we calculated the form of the ef- 
fective inter-bosonic potential when a synthetic magnetic 
field (SMF) is applied to neutral gaseous Bose-Fermi mix- 
tures. We predict that the fermion-mediated effective 
interaction between bosons has a complicated pattern of 
the frequency dependent magnitude. Moreover, the SMF 
renders the inter-bosonic potential oscillatory with sign 
change, thus switching it between repulsive and attrac- 
tive. As a consequence the resonances appear and BF 
mixture that enters the quantum-Hall regime displays 
surprisingly complex dynamics unreachable in conven- 
tional solid state physics. We expect that our theoretical 
results open up the experimental studies [18] of the renor- 
malized interaction energies in stable many-body phases 
with strong correlations and their dynamical properties. 

Restricting our analysis to the lowest energy band of 
a square optical lattice in synthetic magnetic field, the 
Bose-Fermi quantum gaseous mixture can be modeled 
via the following Hamiltonian [19J: 

T-i = ^ {ubi - 1) - ibijb\bj ~ fj-by^^ mi 

where 6|(cJ) and bj{cj) stand for the bosonic (fermionic) 
creation and annihilation operators; Ubi = h\bi{nfi = 
cjci) measures the corresponding boson (fermion) num- 
ber on the site i, [4 > is the on-site repulsion 
and /i(,(/i/) stands for the chemical potential for bosons 
(fermions). The DD interaction between the bosonic 
and non-interacting, spin-polarized (collisions in the s- 
wave channel are forbidden by their statistics) , fermionic 
atoms is denoted by Ubf and depends, on boson to 
fermion mass ratio mb/mf. Here, identifies summa- 
tion over the nearest-neighbor sites. Furthermore tb(tf) 
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sets the kinetic energy scale for bosons (fermions). 

A synthetic magnetic field B = V x A{r) enters 
the Hamiltonian Eq. ([T|) through the Peierls phase fac- 
tor according to tij — > exp /J'' A • di^ , where 

$0 = he/ e is the fiux quantum and e is an elementary 
charge. Thus, the phase shift on each site is determined 
by the vector potential A (r) and can be controlled ex- 
perimentally [7|. The magnetic field is introduced in the 
theory by the density of states (DOS). There are signif- 
icant difficulties in obtaining and analyzing and analyze 
the solutions of the above analytically for every value of 
/. Only a few closed formulas for DOS are accessible 
|20| and consequently not every applied magnetic field 
can be described theoretically. We expanded the set of 
the available analytical solutions obtaining closed formu- 
las also for / = 1/8 and 3/8. This allows us the detailed 
analysis of the dynamical response function which have 
been found to play a crucial role in complex systems. 
The partition function of bosonic and fermionic mixtures 
is written in the form Z = J [VbVbVcDc] e-'^I^'^^l with 
action given by 

S = Sb+Sc.+ / dTH{T) (2) 
Jo 

where St = J2i Iq '^'^^i-^b., and Sc = Z^i /o'^ '^'^Cj^Cj. 
Using the bosonic (fermionic) path integral over the com- 
plex fields depending on the "imaginary time" < r < 
/3 = l/k^T with T being the temperature we can eas- 
ily integrate over the fermionic fields [U] because spins 
are frozen due to influence of the magnetic trap and 
there is no direct interaction between fermions. After 
that, we obtain the partition function in the form Z — 

J [VbVbVcDc] e-'Si.[fc,»b]e-TrlnG,^ rpj^g ^^^^ ^j^g 

point correlation function for noninteracting fermions Gc, 
after exploiting Fourier-Matsubara transform reads: 

Tr In Gc = - ^ A^ (uje) Xk [ivn) A^k (-w^) , (3) 

where cue = 27re/P {ve = tt {2i + 1) //3) with g = 
0, ±1,±2, ... are the Bose (Fermi) -Matsubara frequencies 
respecting periodic (antiperiodic) boundary conditions 
of the bosonic (fermionic) field operator with A^ (tuf) ~ 
bk (wf ) bk (uJe) and 

^fe(^^^)-E 1/, ' p/, . ' (4) 

is the Lindhard function - more commonly called the ran- 
dom phase approximation with npix) being the Fermi 
distribution; t^i"^ is the dispersion relation calculated 
from Harper equation |20| . It correctly predicts a num- 
ber of properties of the collective phenomena in electron 
gas such as plasmons [5T]. To stay in the local regime we 



perform k and k' integration over the first Brillouin zone 
and, in the T — > limit, using an analytic continuation, 
we obtain imaginary part x" i'^) of the local dynamic 
Lindhard function (LDLF). Therefore, the corresponding 
real part x' can be deduced from the Kramers-Kronig 
relation. We noticed that x" ('^) is proportional to the 
absorption spectrum of the medium so it can be directly 
measured. Now, doing the inverse Fourier transform Eq. 
[3]and using gradient expansion, we obtain quadratic form 
of the trace with extracted frequency dependence 

Tr In Gc ^ - f dr [b, (r) h (r)] ' . (5) 

The consequence of the difference in masses of bosons 
and fermions is the fact that the speed of the Bogoliubov 
sound Ub for bosons differs from the first sound f / of the 
ideal Fermi gas. In typical experimental realizations the 
acoustic long-wavelength boson and fermion velocities are 
comparable and both constituents equilibrates similarly. 
The mentioned different mass ratio has far-reaching con- 
sequences, including the possibility of generating the DD 
oscillations [H]. With these concerns in mind we do not 
restrict our calculations to the static limit but consider 
also the local dynamical response function hence the re- 
tardation effects are an important paxt of the physics. 
Moreover, we see jUj from the C4j {nib/mf) dependence 
that even if the fermion mass is very large system has 
a finite value of the interaction strength. On the other 
hand slow bosonic atoms would affect the system provid- 
ing strong inter-species coupling considerably faster. We 
neglected the back action of the bosons on the fermions 
(renormalization of the fermion DD correlation function) 
because of the diluteness of the system and theory limi- 
tations condition Ub > U^fX' i^) |23j . 

Our analysis is carried out in the frequency domain, 
although measurements are sometimes made in the time 
domain and then Fourier transformed to the frequency. 
When we add Eq. ([5| to the bosonic part of the ac- 
tion there is a striking resemblance to the one-component 
Bose-Hubbard action with the original repulsive interac- 
tion replaced now by Ub — )• t/eff = Ub + U^fX' (<^) which is 
the induced, frequency-dependent, effective inter-bosonic 
potential. From the above we see the DD correlations be- 
tween the constituents give rise to additional interaction 
among bosons, which is robust to repulsive or attrac- 
tive nature of the inter-species interaction but not to the 
sign of the LDLF. Before carrying out an analysis of the 
phase diagrams it is a good idea to examine the struc- 
ture of the response function graphically Fig. [T] Struc- 
ture that will usually reflect, at least in part, the physical 
processes present in the real systems. The interactions 
caused by the DD correlations may change its sign even 
if the system has not the phase shift imprinted on it as 
a result of the collective excitations. For values / = 
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Figure 1: (Color online) Real t[x" d-^) (x-y plane) and imag- 
inary tf-^ {iLi) (z-y plane) parts of the local, frequency uj/Stf 
dependent, density-density response function in the synthetic 
magnetic field for several values of / = 0, 1/2, 1/3 and 1/6 (for 
the case with / = 1/4 see supplemental material [24]). The 
complex plot (x-z plane) is the "Cole-Cole" like diagram |25| . 
The three-dimensional parametric curve shows the evolution 
of the density-density response function with the frequency 
and guarantees the stability of the system according to the 
Nyquist theorem (see text). Normalization here tf = t^^''^". 
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Figure 2: (Color online) Real t/x" {^) and imaginary tfx' (<^) 
parts of the local, frequency aj/8t/ dependent, density-density 
response function in the synthetic magnetic field for / = 1/8 
and / — 3/8. Normalization here tf = t^J''^^ . 



and 1 /2 we have only one resonance and induced attrac- 
tive interaction between bosons becomes repulsive above 
some critical frequency Wcrit- Decreasing the strength 
of SMF makes the physics nontrivial and for g > 2 the 
number of the resonances increase significantly Fig. |2] 
Such complex behavior of the mixtures emerges in the 
limit not reachable in conventional systems of condensed 
matter physics because the very high values of magnetic 
field are required to acquire the desired range / < 1/2. 
The method of an optically synthesized magnetic field for 
ultracold neutral atoms |7j already reached that regime 
with a high accuracy. We see that the situation is greatly 
modified if one applies lower magnetic fields. The nar- 
rower peaks in frequency domain the longer time we need 
to observe the collective responses from the system. Ac- 
cording to the above we notice that the BF mixtures has 
to be in a very high magnetic field to make oscillations 
experimentally observable and seen as a density modula- 
tion. The stability of the effects we analyze comes from 
the Nyquist theorem for the frequency dependent dia- 
gram Fig. [T| The parametric curve x' (w) vs x" i'-^) ^ 
a function of frequency and the ensuing effective interac- 
tion described by its behavior t/off automatically yields a 
stability if the curve does not encircle the origin |26| . 

To obtain an equation of state we apply the quantum- 
rotor approach, that successfully casts an essential part 



of the physics of strongly interacting fermions 
and bosons [3D] exclusively, to BF mixtures. We incor- 
porate fully our calculations to the phase fiuctuations 
hi (r) — [ao + a- (t)] e*"^'^"^' governed by the gauge U (1) 
group and drop corrections to the amplitude (t) of 
the order parameter V^b = (5j (r) e*"^'^"^^) = feoV^B- The 
non- vanishing value of the \1/b signals a bosonic conden- 
sation inni- A phase 4)i (r) of the many-body wave func- 
tion might be arbitrary but correlations among the local 
phases of its constituents can bring unusual gauge struc- 
tures. Now, the partition function Z = J [!?(/)] e~'^p'»['^l , 
with an effective action can be expressed in phase-only 
terms 



^ph [(!>] 



dr 



E 



2C/cff iUcS 



4>^ (r) 



<}pi[<t>i{-r)-'l>j{r)] 



(6) 



The phase stiffness coefficient given by Jp/« 
tl^'^ (stl^'' + flb-UbfNF^ /Ub describes the hopping 
matrix elements renormalized by the amplitude of the or- 
der parameter and N-p is the average number of fermions. 
The critical line equation that separates the Mott in- 
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sulator - superfluid transition Fig. [3j details of similar 
derivation of the critical line equation are described in 
|14| . will take simple form: 

where ^ = ^S/ax — C ^^^id ^S/ax is the maximum of the 
band spectrum. The renormalization parameters are de- 
fined as: a = 1 + (w) and ^ = ^iVp - 1/2. In Eq. 
([t]) u iji/U) ~ frac (fi/U) — 1/2, where frac (a;) is the frac- 
tional part of the number. Because the higher values of 
the normalized chemical potential for the fermions Hj/tf 
decreases x' (w) and x" {^) jl4|, consequently terms con- 
taining explicitly the average density of fermions N-p 
will acquire more significance than these with exclusively 
the inter-species interaction Ubf- The periodicity of the 
phase diagram can be easily characterized by its evolu- 
tion with changing the a parameter and applied magnetic 
field (see Fig. [3|. For / = 0, taking u! < Wcrit) we recover 
the previous theoretical results in which, after adding 
fermions to the system, the effective interaction UcS (w) 
becomes smaller than repulsive energy Ub for bosons only 
and superfluid phase increases. However, in the local dy- 
namic limit, when lu > Wcrit the Mott insulator phase 
becomes stronger and bosons tend to localize on the lat- 
tice sites. Imprinting a phase factor on neutral particles 
/ 7^ can provide very complex behavior of the mix- 
tures and induced effective interaction UcS can oscillate 
from attractive to repulsive depending on frequency. The 
measurements by the Bragg spectroscopy method of the 
linear response of correlated two-dimensional BF mix- 
tures at non-zero momentum transfer can give us more 
insight into the excitations spectra modified by the ap- 
plied synthetic magnetic field. 

In conclusion, we have studied a planar mixture of 
bosons and spinless fermions with synthetic magnetic 
field imposed on the system. We found that the un- 
derlying dynamics of mixture of particles with different 
statistics and masses entering a quantum-Hall regime is 
very complex allowing effective bosonic interaction to be 
switched between repulsive and attractive. The experi- 
mental evidence of our findings is feasible however precise 
measurements of the magnetic field are requisite which 
is possible with the recently developed optically synthe- 
sized magnetic field for neutral atoms |7]. 

We thank T. K. Kopec, R. Micnas and I. Spielman for 
discussion and comments regarding the paper. 



[1] I.Coddington et ai, Phys.Rev. A 70, 063607 (2004). 

[2] S.Tung et al, Phys.Rev.Lett. 97, 240402 (2006). 

[3] V.Schweikhard et ai, Phys.Rev.Lett. 99, 030401 (2007). 




Figure 3: (Color online) The phase diagram showing the Mott 
insulator (MI) - superfluid transition (MI phase inside the 
lobes, superfluid above the surface) obtained from Eq. [t] of 
the Bose-Fermi mixtures confined in two-dimensional square 
optical lattice with synthetic magnetic field / = 1/2 in the 
space of the parameters ti,/Ub-i-ib/Ub-a for 77 — —2 {uj > Wcrit 
means that system is in the repulsive regime (7cff > Ub with 
x' {(jj) > 0, even though the interspecies interaction is negative 
Ubf < 0, see also Fig. [T]). Upper panel is the density plot of 
the three-dimensional phase diagram. 



[4] R.Bhat et ai, Phys.Rev. A 76, 043601 (2007). 

[5] N.R.Cooper, Advances in Physics 57, 539 (2008). 

[6] A.Leggett, Quantum liquids (Oxford, New York, 2006). 

[7] Y.-J.Lin et al, Nature, 462, 628 (2009). 

[8] C.Ospelkaus et ai, Phys.Rev.Lett. 96, 020401 (2006). 

[9] F.Ferlaino et al, Phys.Rev. A 73, 040702(R) (2006). 
[10] K.Ciinter et ai, Phys.Rev.Lett. 96, 180402 (2006). 
[11] Th.Best et ai, Phys.Rev.Lett. 102, 030408 (2009). 
[12] J.Catani et al, Phys.Rev. A 77, 011603(R) (2008). 
[13] R.M. Lutchyn et ai, Phys.Rev. A 79, 011606(R) (2009). 
[14] T.P.Polak, T.K.Kopec, Phys.Rev. A 81, 043612 (2010). 
[15] P.Buonsante et al., Phys.Rev.Lett. 100, 240402 (2008). 
[16] G.Mazzarella, Eur.Phys. J. D 50, 61 (2008). 
J17] F.Illuminati, A.Albus, Phys.Rev.Lett. 93, 090406 (2004). 
[18| S.Will et ai, Nature, 465, 197 (2010). 
[19[ A.Albus et al, Phys.Rev. A 68, 023606 (2003). 
[20[ T.P.Polak, T.K.Kopec, Phys.Rev. A 79, 063629 (2009). 
[21[ G.Giuliani, G.Vignale, Quantum Theory of the Electron 

Fluid, (Cambridge, 2005). 
[22[ A.Mering, M.Fleischhauer, Phys.Rev. A 81, 011603(R) 
(2010). 

[23[ D.V.Efremov, L.Viverit, Phys.Rev. B 65, 134519 (2002). 
[24[ See supplemental material at 

|http: / /link.aps.org /supplemental/1 10.1103/Phys- 



5 



RcvA. 83. 011602]. 
[25] K.S.Colc, R.H.Cole, J.Chem.Phys. 9, 341 (1941). 
[26] R H. Chavanis et al, Eur.Phys.,J. B. 69, 389 (2009). 
[27] S. Florens, A. Georges, Phys.Rev. B 70, 035114 (2004). 



[28] T. K. Kopec, Phys.Rev. B 70, 054518 (2004). 

[29] E. Zhao et al, Phys.Rev. B 76, 195101 (2007). 

[30] T.P.Polak, T.K.Kopec, Phys.Rev. B 76, 094503 (2007). 



